Quantum transport through ballistic cavities: soft vs. hard quantum chaos 
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We study transport through a two-dimensional billiard attached to two infinite leads by numerically 
calculating the Landauer conductance and the Wigner time delay. In the generic case of a mixed 
phase space we find a power law distribution of resonance widths and a power law dependence of 
conductance increments apparently reflecting the classical dwell time exponent, in striking difference 
to the case of a fully chaotic phase space. Surprisingly, these power laws appear on energy scales 
below the mean level spacing, in contrast to semiclassical expectations. 
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Advances in the fabrication of semiconductor het- 
erostructures and metal films have made it possible to 
produce two dimensional nanostructures with a very low 
amount of disorder |Q. At low temperatures, scattering 
of the electrons happens mostly at edges of the struc- 
tures with the electrons moving ballistically between col- 
lisions with the boundary. Theoretical and experimental 
investigations have shown that the spectral and transport 
properties of such quantum coherent cavities, commonly 
called "billiards" , depend strongly on the nature of their 
classical dynamics. In particular, integrable and chaotic 
systems were found to behave quite differently [|]||] • 

Generic billiards are neither integrable nor erg odic @, 
but have a mixed phase space with regions of regular as 
well as chaotic dynamics ||. Their dynamics is much 
richer than in either of the extreme cases, as phase space 
has a hierarchical structure at the boundary of regular 
and chaotic motion. In particular, this leads to a trap- 
ping of chaotic trajectories close to regular regions with 
a probability P(t) ~ t^ 13 for t > t , to be trapped longer 
than a time t, with to of the order of a few traversal 
times [||. The exponent j3 > 1 depends on system and 
parameters with typically (3 ~ 1.5 ||. This power-law 
trapping in mixed systems is in contrast to the typical ex- 
ponentially decaying staying probability of fully chaotic 
systems (see Fig. |l|). 

Recently, it was shown semiclassically employing the 
diagonal approximation that the variance of conductance 
increments (for a small dc bias voltage) over small energy 
intervals AE grows as |7|J|] 

Ag 2 (AE) ee ([g(E + AE) - g{E)f) E ~ |A25|", (1) 

for mixed systems if f3 < 2. This is in strong contrast to 
an increase as (AE) 2 in the case of fully chaotic systems 
||. The semiclassical approximation requires AE to be 
larger than the mean level spacing A, corresponding to 



the picture that quantum mechanics can follow the classi- 
cal power law trapping at most until the Heisenberg time 
tij = h/A In the semiclassical approximation the 
graph of g vs. E has the statistical properties of fractional 
Brownian motion with a fractal dimension D = 2 — (3/2 
0. Fractal conductance fluctuations have indeed been 
found in experiments on gold wires |ic| and semiconduc- 
tor nanostructures |ll| and numerically for the quantum 
separatrix map jl2| . 

In this Letter, we numerically study quantum trans- 
port through a simple cavity, the cosine billiard [jl3| (see 
insets of Fig. [j]). Although we observe completely differ- 
ent behavior for the mixed and fully chaotic cases in the 
semiclassical regime of many (45) transmitting modes, we 
find in the mixed case no indication of fractality or frac- 
tional Brownian motion behavior of the graph g vs. E. 
This is the first surprise, as it is in contrast to the above 




FIG. 1. Classical dwell time probability P(t), for the cases 
of mixed (thick line) and fully chaotic (thin line) dynamics 
with t in units of the traversal time. The insets show the shape 
of the billiard with attached leads and a Poincare surface of 
section of the closed billiard for the mixed (upper right) and 
chaotic (lower left) case. 
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mentioned semiclassical JjJ , experimental , and nu- 

merical E^] works. Instead, the conductance is charac- 
terized by narrow isolated resonances, with the classi- 
cal exponent f3 appearing in a power law distribution 
of resonance widths smaller than the mean level spac- 
ing. This leads to a scaling of Ag 2 (AE) in agreement 
with the semiclassically derived Eq. ([!]), however, only 
on scales below the mean level spacing. This surprising 
result contradicts the semiclassical intuition that quan- 
tum mechanics may mimic classical properties at most 
until the Heisenberg time corresponding to energy scales 
above the mean level spacing. At present, there is no 
explanation for these numerical results. They show that 
even with a detailed (semiclassical) knowledge of the uni- 
versal chaotic regime as well as the integrable case at 
hand we are just at the beginning of understanding the 
quantum properties of generic Hamiltonian systems. 

The cosine billiard (l3) is defined by two hard walls 
at y = and y(x) = W + (M/2)(l - cos(2ttx/L)), for 
< x < L, with two semi- infinite perfect leads of width 
W attached to the openings of the billiard at x = and 
x = L (see insets of Fig. [[]). By changing the parame- 
ter ratios W/L and M/L the stability of periodic orbits 
associated with the billiard can be changed, allowing a 
transition from a mixed to a predominantly chaotic phase 
space. Note, that in the mixed case the leads couple to 
the chaotic part of phase space only. 

The ^-matrix of the system has been calculated by the 
recursive Green's function method after expanding the 
two-dimensional wave function in terms of local trans- 
verse energy eigenfunctions [ pT[ . In the numerical calcu- 
lations, it was checked that a sufficient number of modes 
in the expansion in transverse eigenmodes was kept and 
that the lattice constant in x-direction was sufficiently 
small. For a given energy Ef — h 2 k F /2m, N modes 
in the leads are transmitting, with kpW/ir > N. We 
turn from the i'-matrix to the experimentally relevant 
conductance at small dc bias voltage using the Landauer 
formula, G — e 2 /hTr (tt^), where t is the transmission 
matrix. Spectral information is contained in the Wigner- 
Smith time delay r = -ihTr (SUS/dE)/2N, where 2N 
is the dimension of the S'-matrix. All energies in this 
paper are given in units of h 2 n 2 / (2mW 2 ) . 

Figure || shows the dimensionless conductance g = 
G(h/e 2 ) and the Wigner-Smith time delay r [in units of 
2mW 2 1 (Tiit 2 )] for parameters corresponding to a mixed 
phase space (W/L = 0.18, M/L = 0.11) and a chaotic 
phase space with no apparent stability island (W/L — 
0.36, M/L = 0.22) for N = 45 transmitting modes. The 
differences are quite dramatic. For the fully chaotic case, 
both quantities are smooth functions of energy and in 
good agreement with semiclassical theory (see below). 
While the average values are comparable, many sharp iso- 
lated resonances on top of a smooth background are visi- 
ble in the mixed case 15 1 , also in contrast to the semiclas- 
sically predicted fractional Brownian motion. The simple 
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FIG. 2. Dimensionless conductance g (a) and Wigner time 
delay r (b) for N = 45 propagating modes. Thick (thin) lines 
for mixed (fully chaotic) case. 



explanation that these narrow resonances are related to 
quantum tunneling into the islands of regular motion |l(j| ] 
does not apply here, as the phase space volume of stable 
islands is about 5%, while the narrow resonances (below 
the mean level spacing) make up about 18% of all states 
associated with the billiard. This roughly corresponds to 
the phase space volume around the stable islands where 
trapping of chaotic trajectories occurs. 

In order to analyze the narrow resonances in the mixed 
case, it is convenient to examine the Wigner-Smith time 
delay. Each resonance in the time delay has the Breit- 
Wigner shape, characterized by a width Ti and a height Tj 
situated at an energy Ei on top of a smooth background. 
We find our data well described by 



r(E)=J2- 



rf/4 



(E-Etf+Tt/1 



^"smooth (E), (2) 



with T smoot h(-£') oc E~ x l 2 . Since the phase shift through 
a resonance is 27T, width and height are related by r^Ti = 
2/N. The energy was initially sampled on an equidis- 
tant grid and subsequently refined in order to resolve the 
sharp resonances. Only resonances with a V < 10 -3 , i.e. 
much smaller than the initial grid are lost. As a result 
we can numerically construct the cumulative distribution 
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FIG. 3. Normalized cumulative distribution of resonances 
widths N(T) for the mixed case. The thin line serves as a 
guide to the eye. The mean level spacing A is shown in the 
figure. 



N(T) of resonance widths, corresponding to the proba- 
bility of finding a resonance smaller than T (Fig. ||). The 
distribution is very broad, spanning 5 orders of magni- 
tude, and is approximately a power law N(T) rj aT r , 
with r ~ 0.35, over a wide range below the mean level 
spacing A = 0.176. 

The consequences of this broad distribution of reso- 
nance widths for the variances of conductance and time 
delay increments are studied now. For AE <C A correla- 
tions between different isolated resonances (r <C A) do 
not contribute to the variance Ag 2 , which then is gov- 
erned by the distribution of resonance widths. Each res- 
onance is reflected in the conductance, 

n r 

g{E) = 5smoot h(£) + 5 9i( E )> (3) 

i=l 

where Sgi (E) is a function of the width Ti and the typical 
height Nr is the number of resonances with < 
A in the energy interval Eb — Ea over which we take 
the average. The variance of the increments of a single 
resonance is given by 

({8g l (E + AE)-5g l (E)} 2 ) E = 

g|r t f hiAE/Ttf , ae « r 4 

E B - E A \ 1 , AE » T 4 ' [V 

which defines <?i and where bi is a numerical factor of or- 
der unity. Since the distribution of widths is very broad, 
the strong inequalities are almost always fulfilled in the 
sum over resonances. Splitting this sum into resonances 
smaller and larger than AE, we get 

\r 4 <AE Ti>AE 1 / 

(5) 



Replacing the sums by integrals over the density of 
widths n(r) w arT r ^ 1 and neglecting the weak fluctu- 
ations of gi and bi as compared to Tj, we can estimate 
for small AE, 

Ag 2 (AE) ex G 2 ) EB N * EA a\AE\ 1+r . (6) 

where (• • •) stands for the average over isolated reso- 
nances. A power law distribution of resonances thus leads 
to a power law increase of the variance of conductance 
increments with the exponent given by 1 + r. 

Fig. [| a) shows the variances of the conductance incre- 
ments. On scales smaller than the minimum resonance 
width the variance increases quadratically, as expected. 
On larger scales we find the power law Eq. (||). This 
result coincides with the semiclassically derived Eq. (|l|) 
with r — (3 — 1, however, only on scales below the mean 
level spacing. At present, there is no explanation why 
the classical exponent (3 appears on such small energy 
scales. Remarkably on scales above the mean level spac- 
ing the correlation energy for the conductance fluctua- 
tions is given by the Weisskopf width Tyy » 2, as in the 
fully chaotic case (see below). 

The variance of increments of the time delay are shown 
in Fig. [| b). Since the variance At 2 (AE) measures the 
square of the resonance peak height in the time delay, in 
the mixed case, they are completely dominated by the 
sharpest resonance, once the energy exceeds the mini- 
mum resonance width. Thus, in contrast to fully chaotic 
systems, in mixed systems the scale of the correlations of 
the time delay is the smallest resonance width and not 
the Weisskopf width Yw- 

For comparison in Fig. ^ we also show the results 
for Ag 2 {AE) and At 2 {AE) for the fully chaotic case. 
They are characterized by single scales T g m 4.8 and 
r r w 3.8 and are in good agreement with semiclassical re- 
sults ||[l8| . The chaotic case can also be described by the 
random matrix theory (RMT) (T^| , whose results coincide 
with the cited semiclassical ones for N 3> 1 [^(|. In the 
absence of direct processes, random matrix theory pre- 
dicts a single correlation scale, known as Weisskopf cor- 
relation width, Tw = A/27r53 c T c , where the sum runs 
over all channels c with transmission probability T c [ ^l) . 
Approximating ^ c T c by twice the average dimensionless 
conductance we obtain ~ 4.2, in agreement with the 
numerical values within the statistical accuracy. Before 
concluding, it is worthwhile to stress that depending on 
N and the coupling to the leads, quantum chaotic scat- 
tering can also exhibit isolated resonances. Their width 
distribution, however, follows a ^-distribution with 
degrees of freedom M], rather than power law. 

In conclusion, we have shown that generic Hamilto- 
nian systems, which have regular as well as chaotic phase 
space regions, differ drastically in the Landauer conduc- 
tance and Wigner time delay from fully chaotic systems. 
We find many isolated narrow resonances with a power 
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FIG. 4. Variance of the increments of the conductance 
Ag 2 (AE) (a) and Wigner time delay Ar 2 (A_E) (b). The thick 
(thin) lines correspond to the mixed (fully chaotic) case. The 
thin dashed lines are fits to the semiclassical expressions for 
the chaotic case. The thick dashed line in Fig. b) is the con- 
tribution of the sharpest resonance to the variance for the 
mixed case. For comparison, the ticks on the upper borders 
mark the widths of the individual resonances. The inset in a) 
shows the corresponding data for a quantum graph modeling 
a mixed phase space (f3 = 1.48) and was provided by the au- 
thors of ref. jlTj. The mean level spacing A is shown in the 
figures. 



law distribution of their widths accompanied by a power 
law increase of the variance of conductance increments. 
Both power laws appear to be connected to the classi- 
cal power law trapping, surprisingly they only appear on 
scales below the mean level spacing. Similar unexplained 
power laws are found in recent studies using quantum 
graphs p2| m odeling a mixed phase space (see inset of 
Fig. ^ a) ]17fl. Further research on the quantum signa- 
tures of classically mixed systems is urgently needed. 
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Note added. — The authors of Ref. have informed 



us that there are quantum graphs where the classical and 
quantum exponents do not agree. 
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